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^ ABSTRACT 

> : I„ order to flud a law characterising the decrease of velocity along a jet. five 

O ■ analytical methods are suggested. The first two simple models examine the variation 



of velocity in the presence of Newton's or Stoke's resistance. The equation that 



^ I represents the conservation of the momentum along a pyramidal sector is solved from 

• an analytical point of view (third model). The application of the conservation of the 

Q ■ total momentum flux allows us to deduce the velocity of the galaxy as a function of 

H ; 

c/3 ! time for classical velocities (fourth model) and relativistic velocities (fifth model). 
■ 

The variation of velocity along the jet combined with an adequate composition of jet 



^ ■ precession velocity , rotational velocity of the galaxy , and galaxy dispersion velocity 

in the cluster allows us to trace the geometrical pattern of the head-tail radio 
sources. Application of the developed theory/code to the radio galaxies 
NGC1265,NGC4061,NGC326, and Cygnus A gives the central galaxy's approximate 
dispersion velocities in the direction perpendicular to the jet. A transition from 
head-tails to classical double radio galaxies as a function of the increasing jet's 
mechanical power is introduced. 
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1. INTRODUCTION 



The extra-galactic radio sources are classified on the basis of the position of the 
brightest radio emitting regions with respect to the channel: FR-I have hot spots 
that are more distant from the nucleus (typical example Cygnus A) , FR-II radio 
galaxies have emission uniformly distributed along the channel (typical example 
3C449). The physical parameter that governs this classificatio n is the radiated p ower: 



Fanaroff fc Riley 



the na ore powerful radio galaxies being classified as FR II,see 
( 1974^ for details. The classification then becomes more complicated in the presence 
of complex morphologies and two classifications are introduced : NAT ( Narrow 
Angle Tail) and WAT ( Wide Angle Tail). But a closer look at NGC 1265 reveals 
that the angle subten ded by the two channels is not so " narrow" , see Figure 1 of 
Odea fc OwenI (|l986| ): possibly, in this case a distinction should be made between 
internal regions of the channel where the patterns are similar and the extended 
regions where conversely the division in WAT and NAT is more evident. 
For more details on the classification scheme the interested reader is advised to refer 
to the book by 



de Yound (|2002i ). The X-ray observations of the Galactic black hole 
XTE Jl 118+480 are consistent with extended jets being the source of the hard X-ray 
flux. In this hypothesis the disc w ould then simply represent a small so l id aii gle as 



seen from the emission region, see 



Miller. Ballantvne. Fabian, fc LewinI (120021). T his 



point of view is not new and has been explored bv lMendoza. Hernandez, fc Lee 
(20051) who suggested a unified model for quasars and mu-quasars . In general, the 
shapes seem to be independent from (or weakly correlated to) the cluster 
parameters, such as physical position in the cluster, number of galaxies, cluster 
richness, cluster morphology and x-ray emission. Some of the models that can 
explain the bending of the jets will now be briefly reviewed 



The geometrical models started with 



Jaffe fc Perolal (|l973( ) where two models 



were developed to explain the formation of "tailed" radio sources like 3C 129. 
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Th ey continued inserting 



see 



p recess ion and relativistic effects in SS433, 



Hiellming fc JohnstonI (|1982[ ). An analytical model was prese nted for th e 



evolut ion of a powerful double radio source on a small scale , see 



(200' 



Alexander 



Euler's equation model, where bending is produced b y the ram pressure 
PicuVj of the intra-cluster medium, is analysed by iBurns fc OwenI (|l980( ) in 



order to explain the structure of 1638+538 (4C 53.57), and by 



Venkatesan. Batuski. Hanisch. fc Burm 



( 1994( 1 in order to understand the 



behaviour of NAT in poor clusters of galaxies. Further on, a hydro-dynamical 

code was built in order to explain the NAT sources, see 

Norman. Balsara. &: O'Deal (|l994^ . From an observational point of view it is 

interesting to note that the velocities of the NAT galaxies are inadequate for 

producing the ram pressure necessary to ben d the radio jets, see 

Bliton. Rizza. Burns. Owen. &: Ledlowl (|l998h . 



The slingshot ejection model consists in the ejection of black holes from the 
host galaxy: the bending is obtaine d from the oscillations about t he center of 



the merged galaxy, see for example 



Valtonen &: KotilainenI ( 198' 



The 3-dimensional magneto-hydr odynamic (MHD) simulations 



3ased on the 



Sweeping Magnetic-Twist model 



Nakamura. Uchida. fc Hirosd (|2001l ). produce 



a wiggled structure of AGN radio jets. 



Raga fc Cantol mm 



Canto &: Raga 



The interaction between a jet with a stratified cloud, see 
or wit h a spherically symmetric pressure stratification , see 
( 19961 ). can lead to a final configuration in which the jet has bored a hole 
through the cloud, or large deflections are obtained. 



A jet can generate an internal shock wave by adopting either a non-relativistic 
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framework , see 


Icke & Huehes 


1 1 

Mendoza & Lonsraiii ((2002 


)• 



From a relativistic point of view 



Mendoza fc Longaiii (j2001l ) showed that 



bending relativistic jets is much more difficult than bending non-relativistic 
ones. This is the reason why most FR-II radio galaxies appear straight. 

The already cited models concerning the bending leave a series of questions 
unanswered or partially answered: 

• Is it possible to deduce a law of motion in the presence of viscosity in laminar 
and turbulent jets? 

• Is it possible to trace the jet velocity on the basis of physical parameters such 
as linear density of energy of the radio source and constant density of the IGM 
(Inter Galactic Medium)? 

• Is it possible to develop a consistent theory in which at a given time the jet 
opening angle increases abruptly? 

• Is it possible t o include relativistic e ffect s in turbulent jets tha t generalise the 



work made by 



Raga fc Cantol (|199a) and 



Canto fc Raga 



( 19961) on the 



interaction between a jet and a stratified cloud ? 

Once the main jet trajectory is obtained, could other kinematical components 
be added such as jet precession, galaxy rotation and galaxy bulk velocity? 

Does the developed theory match the observed patterns traced by the radio 
sources? 



Could the transition FR-I =^ FR-II be simulated by increasing the total jet's 
luminosity? 
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In order to answer these questions in Section |21 the equations of motion were derived 
in the presence of Stoke's law of resistance , Newton's law of resistance and turbulent 
eddy viscosity. Some analytical computations on the momentum conservation in an 
extra-galactic pyramidal sector advancing in a surrounding medium were developed 
in Section ini In Section H] the theory of the two-phase beam was set up. The theory 
of the relativistic turbulent jets was developed in Sectional The theory for the 
composition of the velocities in a head-tail radio source was developed in Section ^ 
The application of the developed theory to well specified radio sources such as 
NGC1265,NGC4061, NGC326 and Cygnus A was reported in Section □ A power 
transition is simulated in Section 1731 replacing the energy with the total power in the 
law of motion. The theory of the Kelvin-Helmholtz instabilities was reviewed and 
implemented on the radio-sources in Section El 



When a laminar jet moves through the IGM or the ISM a retarding drag force Fdrag 
is applied. If v is the instantaneous velocity the simplest model which is usually 
considered is where 



with n as an integer. Here the case of = 1 and n = 2 is considered. In classical 
mechanics n = 1 is referred to as Stoke's law of resistance and = 2 as Newton's law 
of resistance. We now consider two laminar jets and the turbulent case separately. 



2. VISCOSITY MODELS 



n 




2.1. Stake's behaviour 



Consider a slice of the jet with mass m. The equation of motion is given by 



dv 

m— = — 



(2) 



dt 
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or 

I— ' 

where B = b/m. The temporal behaviour of the velocity turns out to be 

v = Voe-''' , (4) 

where Vq is the initial velocity. The distance at the time t is 

x = |(l-e-^*) , (5) 

and this can be considered the first law of motion. The velocity space dependence is 

V — Vq = —Bx . (6) 

The velocity of the jet is allowed to vary between a maximum value at the beginning 
(the inner region), v = Vq = f-^ with c/ a number bigger than one , and a minimum 
value at the end (the outer region) , for example v = yqUoo' ^^^^ coefficient 
B can be found by identifying x with the jet's length xi 

^ ^ Vq - c/10000 

Xi 

The lifetime of the radiosource ■,t^s^'^i can be derived from equation (j^ once the 
numerical value of B is known 

/ 1)0X10000 

tTs'' = ^^^ , (8) 

that becomes 

10000 

4^ = ^ • (9) 
2.2. Newton's behaviour 

The equation of motion is now 

m— = —av , (1(J) 
at 
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or 

where A = a/m . The temporal behaviour of the velocity turns out to be 

- = T^ ' (12) 

and the distance 

. = , (13) 

which can be considered the second law of motion. The velocity space dependence is 

x = — ^ • (14) 

The coefficient A can be found by inserting ^ = 1/10000 which corresponds to 

f = 30 km/sec while L, the length of the jets, is equal to x , later called xi. 

Here the lifetime of the radiosource ,^^1*"*°" , can be derived from equation (jl2p once 

the numerical value of A is known 

2.2.1. Parameters of the artificial viscosity from a sample of WAT 



A sample of 7 radio-galaxies clas sified as Wide Angle Ta i 



sensitively at high resolution by 



s (W AT) was imaged 



Hardcastle fc Sakelliod (j200J). This sample is here 



considered as a source of data that allows us to derive the parameters of Section 127 
and Section 12. H see Table ^ The parameters of the first two laws of motion are 
found once the averaged velocity at the beginning of the jet (alias Cf ) and at the 
end of the jet are fixed by observational arguments. The spectroscopical 
measurements are a good candidate to find the jet's velocity but attention should be 
paid to the fact that at a given position x the velocity can vary considerably from 
the boundary to the center of the jet , see formula (j^D). 



8 



ZANINETTI 



TABLE 1 

Artificial viscosity parameters of seven WAT ; 
Cj — 1.25 and velocity at the end v/c=l/10000. The value of cj is chosen in order to 
have relativistic velocities in the initial stage (inner region) of the jet . 



Radio name xi [kpc] tf^'^'^flOVar] S[l/year] t^|'"*°'^[10Vear] ^[1/pc] 



0647 + 693 


81 


2.9 


3.10 10-6 


2.810^ 


1.1 10- 


1231 + 674 


35 


1.25 


7.13 10-6 


1.210^ 


2.56 10" 


1333 + 412 


20 


0.79 


1.24 10-5 


71 


4.49 10- 


1433 + 674 


49 


1.76 


5 10-6 


1.710^ 


1.83 10" 


2236 - 176 


44 


1.58 


5.67 10-6 


1.510^ 


2.0 10- 


3^465 


28 


1.0 


8.9 10-6 


99 


3.2 10- 


1610 - 608 


13 


0.46 


1.8 10-5 


46 


6.9 10- 
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2.3. Turbulent jets 

The theory of turbulent jets emer ging from a circular hole can be found in different 



Bird. Stewart, T.ightfootI (j2nn2[). Il.a^ndaul (jl987l) 



boo ks with different theories , see 

and iGoldsteinI (|l965l ) . The basic assumptions common to the three already cited 



approaches are 

1. The rate of momentum flow , J, represented by 



J = constant x pb'^vl , 



(16) 



is constant; here x is the distance from the initial circular hole , b{x) is the jet's 
diameter at distance x, Vx,max is the maximum velocity along the the 
centerline, constant is 



constant = 27r / f 
Jo 



where 



with ^ 



X 



h 

'^x,max "l/2 

and p is the dens ity of the surrounding medium see equation (5.6-3) in 



Bird et al. 



( 20021) 



(17) 



2. The jet's density p is constant over the expansion and equal to that of the 
surrounding medium. 

When omitting the details of the computati on, an expression can be found for the 
average velocity Vx-, see equation (5.6-21) in 



BiT'fl f^t al 



mm. 



V 



(t) 



Vx 



2Ci 



^ [i + i(c^3^)^]' ' 

where z/^*-* i s the kinematical eddy viscosity and C3 , see equation (5.6-23) in 



(19) 



Bird et al. 



( 20021 ) 
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An important quantity is the radial position , r = 61/2 , correspo nding to an axial 
velocity one-half of the centerline value , see equation (5.6-24) in 

Vx{hl/2,X) _ 1 _ 1 



Bird et al. 



[X) 



[i + i(^3^ 



(2002), 

(21) 



The exper i ments i n the range of Re ynol ds number ,Re, 10^ < Re < 310^ (see 



Reich ardt 



(IMa) ■ iReichardtl (jlPSlh and 

5i/2 = 0.0848a; 



SchlichtingI ()l979l )) indicate that 



and as a consequence 



and therefore 



Co = 15.17 



(22) 



(23) 



1 



[1 + 0-414(^)2] 



'Vx,max y ) 

Figure 2] reports the behaviour of the velocity distribution as given by equation 
The average velocity, Vx, is ~ 1/100 of the centerline value when r/61/2 = 4.6 and 
this allows to say that the diameter of the jet is 



(24) 



6 = 2 X 4.661/2 • 

On introducing the opening angle a , the following new relationship is found 

a 4.661/2 

— = arctan — 

2 X 



(25) 



(26) 



The generally accepte d relationship be tween the opening angle and Mach number 
see equation (A33) in lde Yound (|2002| ) , is 



a Cs i- 
2 ~ Vj ~ M 



(27) 



where Cg is the sound velocity, Vj the jet's velocity and M the Mach number. The 
new relationship replaces the traditional relationship ()27|) . The parameter 61/2 
can therefore be connected with the jet's geometry 



7 1 /(^\ 

61/2 = — tani— )x . 
/ 4.6 ^2' 



i2t 




Fig. 1. Mean velocity profile vs. channel radius in a circular jet in turbulent flow. The ve- 
locity distribution is a function growing from zero (at the wall of the channel) to a maximum 
value in the central region. 

If this approximate theory is accepted, equation (f^^ gives a = 42.61° : this is the 
theoretical value that originates the so called Reichardt profiles. The value of 61/2 
fixes the value of C3 and therefore the eddy viscosity is 



13 J 1 / 3 1 

i^*-*^ = \ T7^\ -77- = \ TT^V constant hv^^rnax-pr ■ (29) 
V loTT \ pC^ V lovr 63 

In order to continue the integral that appears in constant should be evaluated, see 
equation (fTTj) . The numerical integration gives 

/ /2erfe = 0.402 , (30) 
Jo 

and therefore 

constant = 2.528 . (31) 

On introducing typical parameters of jets like a=5°, Vx,max=vioo = f [km/sec]/100, 
b = bi where 61 is the momentary radius in kpc , it is possible to deduce an 
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astrophysical formula for the kinematical eddy viscosity 

= 0.0028 Mioo^J^ when C3 = 135.61 . (32) 
10 'year 

This paragraph concludes underlining the fact that in extragalactic sources it is 
possible to observe both a small opening angle , ~ 5°, s ee Section 17.41 and gre at 
opening angles , i.e. ~ 34° in the outer regions of 3C31 



Laing fc Bridle! (|200J). 



3. A CONTINUOUS MODEL 

The jet is now explained as due to a continuous flow in a given direction. 
The constrained basic elements of our theory are the mean spread rate of the jet 
measured on the radio maps , the linear density of energy and the initial jet radius 
Tj. The turbulent hyper-sonic flow , Mach number > 6 , is one of the most 
complicated problems in fluid mechanics. Here the key assumption is that due to the 
turbulent mass transfer, the density in the iet is the same as the surroun ding 



medium, see for example equation (10.25) in iHughes fc BrightonI (|1967^ . This 
allows the jet law of motion to be expressed in terms of the IGM's density. We 
identify our jet with a pyramidal section characterised by a solid angle A VL and 
overall length Xi. From a practical point of view, (the first opening angle) is 
reported in the captions. The solid angle dVt in spherical coordinates (r,^, 0) is 

dVL = sm{e)ded(j) . (33) 

On shifting to the finite differences, the following is obtained 

27r 27r 

A^] = (cos(O) - cos(ao^))(tto^) steradians , (34) 

3o0 3o0 

that when the angles are small becomes 

AQ ^ -(— — )^steradians , (35) 
2^360°^ ' ^ ^ 

where ao represents the limit of the integration and is expressed in degrees. 
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Fig. 2. Two dimensional view of the continuous jet; the physical parameters are the same 
as the simulation of NGC4061 , see Figure El 

As an example, when ao=10°, A Q = 2.65 10^'^ steradians. A plot that depicts the 
assumed geometry is reported in Figure |21 

The approximate advancing area of the jet of momentary length x, Aj, will therefore 
be 

Aj ^ x^AQ . (36) 

The jet is divided in many slices perpendicular to the direction of motion , each one 
of thickness (the initial radius); the number of slices, nj, is therefore 

n, = NINT(—) , (37) 

ri 

where NINT denotes the nearest integer. The volume of slice, Vg, is 

Vs = x^Ann , (38) 
and the mass acquired in a slice from the external medium of a density, p, is 

AM = pxVs= px^AQvi . (39) 
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The hypothesis of a jet in which p is constant, is used in the theory of turbulent jets, 
see Section 12.31 The equation that represents the momentum conservation is: 

^ (AMx) = px^Afi , (40) 

where p is the pressure. In the case of constant density, equation becomes 

The pressure of internal gas decreases according to the adiabatic law 

P = P^{^?" , (42) 

where Vi is the volume of the first slice and 

2- 

Under the adiabatic hypothesis, the differential equation (j^T|) is transformed in 

d_ 
di 

where 



(x^x) = i/x"^/^ , (44) 



and 

p, = ^ . 46 

To integrate this equation, x^x = Ax" is used. After adopting the initial condition of 
X = at t = and assuming p is constant, the equation representing the third law of 
motion is obtained 



■l3(3« + 4) . 



3/13 



and since a = 5/6 : 

X = , 
V 90 



/507 \3/i3 
^=(^HtA . (4^ 
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X 



6 X 



The velocity of the jet turns out to be: 

v{t) = 

The astrophysical quantities can now be introduced 



(49) 



3/13 



x{t) = 2.541 kpc 



where E^q =E/10^^ ergs , rj = rj/lkpc 



X, 



x^ 



no 



(50) 



Xi/lkpc , t-i = t/(10'') year and p 



hq mil with no representing the number of particles in a cubic centimetre and mii 
the hydrogen mass. A similar formula is deduced for the velocity 

3/13 



R 



km 

114.92 — 

s 



E56 \ r.! 



i. 



-7/13 



In the case of SNR a comparison can be made with the radius R ( t) oc t^^^ of the 



(51) 



McCrav fc Layzed ( 



Sedovl (|l959|) and 



98 7h ■ and wit h 



LandaiJ (|l987l ): 



adiabatic phase of the SNR , see for example 
R{t) oc t^/^ of the so called Sedov solution , see 

flW = (f — 

4 IT p 

where t represents the time, E is the energy injected in the explosion, and p is the 



(52) 



density of matter. W. 
solution presented in 



len the iets are ana- 



Kaiser fc Alexander 



yzed a comparison can be made with the 
( 19971 ). In their model the density p^ of 



the gas surrounding the jets scales as 



Px = Pol 



d . 



(53) 



where d is the radial distance from the core of the source , po is the initial density , 
Oo is the scale length and f3 a parameter comprised between and 2. Their length of 
the jet Lj scales as 

L.oct^/'-^ , (54) 



and explains the division between FRI and FRII objects in jet power. 
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4. TWO-PHASE CONTINUOUS MODEL 
The startin g point is the c onservation of the momentum's flux in a "turbulent jet" as 



outhned in iLandaul (|1987[ ) (pag. 147). The initial point is characterised by the 
following section : 

AQ = 'Krl . (55) 

Once ao the flrst opening angle , Xq the initial position on the x-axis and Vq the 
initial velocity are introduced, section A at position x is 

A(x) = 7r(ro + (x- Xo)tan((l/2)ao))^ • (56) 

The conservation of the total momentum flux states that 

pvlAo = pv{xfA{x) , (57) 

where v{x) is the velocity at position x . The previous equation is valid if the 
pressure along the jet is constant and viscosity effects are either not considered or 
have the same effect in the whole momentum flux along the jet. Due to the turbulent 
transfer, the density p is the same on both the two sides of equation ()57|1 . The 
trajectory of the jet as a function of the time is easily deduced from equation 



-xq tan ((1/2) ao) + ro - Jrg (r^ + 2 tan ((1/2) Oq) Vq t) 

X = 77— ^ , (58) 

tan ((1/2) ao) 

and this can be considered the fourth law of motion. 

The velocity turns out to be 

v{t) = . (59) 

^ro (ro + 2 tan {{1/ 2) ao)vot) 

In the applications of equation and , r , xq and x will be expressed in 

kpc , the time t in units of 10'' year , vq and v{x) in j^^^^^- The velocity fo can be 
parametrised as a function of the light velocity c as f o = ^ where c/ is an input 
parameter allowed to vary between 1.2 and 10. 
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It is now possible to introduce the two phase beam in which , due to a certain 
physical phenomena, for example the evolution of the K-H instabilities , the beam 
abruptly increases the opening angle that from (first opening angle) becomes ai 
(second opening angle). This phenomena happens at a given time ti ( ti < ^2 ) and 
a corresponding length xi { Xi < X2): here ^2 denotes the age of the radio-source 
and X2 its global length. In this second phase the trajectory is 

-xi tan ((1/2) ai) + rl - ^ rl {rl + 2 tan ((1/2) ai) vi {t - ti)) 

tan ((1/2) ai) ' 

and the velocity 

vit)= ^ . (61) 

y^r, (r, +2 tan((l/2)«i)^; {t-h)) 

5. TWO-PHASE RELATIVISTIC MODEL 
The starting point is the energy momentum tensor ,T'^^, 

T'' = wu'u^ - pg''' , (62) 
where is the 4-velocity , i varies from to 3 , is the enth alpy for unit yo lume , p 



Landadjml). The 



is the pressure and (7*^ the inverse metric of the manifold, see 
condition for momentum conservation in the presence of velocity , v, along one 
direction states that 

iwiy^-^+p)A = cost , (63) 

where A is the considered area in the direction perpendicular to the motion. The 
enthalpy for unit volume is 

w = c^p + p , (64) 

where p is the density , and c the light velocity. The reader may be puzzled by the 
factor in equation (jHSI), where = — However it should be remembered that 
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w is not an enthalpy, but an enthalpy per unit volume: the extra T factor arises from 
"length contraction" in the direction of motion. 

With the assumption of turbulent jets ,p = , the momentum conservation law is 
obtained 

{pv^-^)A = cost . (65) 

J- ,.9. 



This equation is equal to equation (A. 32) , the jet "thrust", in 



de Youne 



Note t he similarity between the previous formula and condition (135.2) in 



20021). 



fll987l ) concerning the shock waves : when A=l they are equals. 
In two sections of the jet we have : 



-J— ^ = p t,2 (^(vr ro^ + vr rg xao + O (oq^) )) ^ ^2 , (66) 



where v is the velocity at position x, Vq the velocity at x=0 , c the light velocity and 
tto the opening angle of the jet. Also here due to the turbulent transfer, the density p 
is the same on both sides of equation (jMj) and the following second degree equation 
in /3 = ^ is obtained: 

+ - p^po^xao - Po^rg = (67) 

where j3o = —. The positive solution is : 



P= I , . (68) 

Y To + xao - i^o a;ao 

From equation (jUHj) it is possible to deduce the distance Xi/2 after which the velocity 
is half of the initial value: 

xi/2 = -3 -1-^ -Y . (69) 



In FR-II's radio-galaxies the velocity is thought to be relativistic until termination at 
hot-spots; formula (j69p states that the typical distance over which the jet is 
relativistic is a function of the three parameters ao , ''"0 and /3o- 
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The trajectory of the relativistic jet as a function of the time can be deduced from 
equation ( and is 

^ -dx = ct . (70) 



^/ ro+xao-Po'^xao 

On integrating it is possible to obtain the equation of the trajectory 

2/3 y-^ ^ ^ ct = . (71) 



ao [-1 + Po ) y/r^Po 
After some manipulation equation (f7T|) becomes 

(-3 ao^po^ + 3 ao^/?/ - Po^ao^ + "o^) + 
(3 ro ao^ -Qrg ao^Po'^ + 3 Po'^ao^) + 
(-3 ro'^Po^ao + 3 ro^ao) x - 
9/4 cH'^ao^ro Po^ - 3 ro^ctao /?o + 3 rg^ctaQ Po^ - 

9/AcH'^ao^ropo'^ + 9/2 cVaoVo = . (72) 

The parameter p , see Appendix |X1 that regulates the solutions of the cubic equation 
is now evaluated. When the equation (f72|) is considered we have p = and therefore 
we have one real root that is (the fifth law of motion ) : 

2/3 



, -2ro + </2^(2ro + ?,ctaoPo-^ctaopo\ , ^ 

x{t) = -\l2 I . o2^ — ■ (73) 



The velocity of the relativistic jet as function of the time is 



rAO - , '''' 3 . (74) 

'2 ro + 3 ctao /3o — 3 ctao Po 



The velocity fo can be parametrised as a function of the parameter Pq , a parameter 
allowed to vary between 0.1 and 0.99999. The two previous equations can be 
expressed in astrophysical units once r\ = ro/lkpc and = t/(10^year) are 
introduced 





{2r^ +9180.0^7^/3 


•q- 9180.0 a po^) 


2/3 


a 1 


[-I+P0') 


1 



x{t) = l/2 ^ ^kpc , (75) 
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Po \R \f2 km 
= 299788.2 , V o v Km ^ ^^^^ 

" 2 +9180.0 tytt/^o - 9180.0^7 a /5o^ ^ 



and 



' Q -^ y±ou.u u: jjQ ^ yiou.u « /j^ 
The transition from relativistic to classical velocity is still represented by 
formulas (jUUI) and but Xi,ri and Vi are the parameters at the end of the 
relativistic phase. 



6. THE CHANGE OF FRAMEWORK 



The wide spectrum of observed morphologies that characterises the head-tail radio 
galaxies could be due to the kinematical effects as given by the composition of the 
velocities in different kinematical frameworks such as decreasing jet velocity , jet 
precession, rotation, and proper v elocity of t he ho st galaxy in the cluster. These 



effects were partially analysed in 



Zaninetti 



()l989l ) ; here part of the developed 



theory and symbols will be used, and will now be briefly reviewed in order to 
represent the law of motion through a matrix. Of particular interest is the evaluation 
of various matrices that will enable us to cause transformation from the inertial 
coordinate system of the jet to the coordinate system in which the host galaxy is 
moving in space. The various coordinate systems will be :K={x,y, z) , 
x^-^^ = {x^^\y^^\ z^-^^) , ... x^^^ = {x^^\y^^\ z^^^). The vector representing the motion of 
the jet will be represented by the following 1x3 matrix 



G 






m 



(77) 



where the jet motion L(t) is considered along the z-axis. 

The jet axis, z, is inclined at an angle '^prec relative to an axis z^^^ and therefore the 
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3x3 matrix, representing a rotation through the x axis, is given by: 



cos (\[^ 



precj 



sin (\[^ 



prec J 



- sin (^prec) COS (^prec) 

From a practical point of view \I/prec can be derived by measuring the half opening 
angle of the maximum of the sinusoidal oscillations (named wiggles) that 
characterises the jet. 

If the jet is undergoing precession around the z^^^ axis, fiprec can be the angular 
velocity of precession expressed in radians per unit time ; flprec is computed from the 
radio maps by measuring the number of sinusoidal oscillations that characterise the 
jet. The transformation from the coordinates x*^^) fixed in the frame of the precessing 
jet to the non-precessing coordinate x*^^) is represented by the 3x3 matrix 



E 



cos {flprec t) siu {flprec t) 
- sin {Qprec t) COS {Qprec t) 











1 



(79) 



The arbitrary orientation of the precession axis relative to the central galaxy at rest, 
necessitates a transformation from the coordinate frame x^^^ to the frame x^^^ . The 
relative orientations are assumed to be characterised by the Euler angles 
(<l>j, Qj, \l'j). There is not uniform agreement on the designation of the Euler angles 
and the manner in w hich they are generated. W e have chosen here to use the 
conventions found in iGoldstein. Poole, fc Safkd (|2002| ) and the 3x3 matrix , De, is 



cos'^jcos^j — cosQjSin^jSin'ifj 
—sin"^ jcos^j — cosQjSin^jCos'^j 



cos'^jsin^j + cosQjCos^jSin'^j sin^> jSinQj' 
—sin'^jsin^j + cosQjCos^jCos'^j cos^ jSinQj 

cosQj 



-sinQjCos^j 
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On assuming that $,= 0° , \l/,= 0° and 9,= 90° we have the simple expression 



D, 



1 



1 



0-10 



Another 3x3 matrix is introduced which represents the transformation from (x^) to 
(x'^) and defined by a rotation through the axis z^^^ of an angle fi^t where is the 
angular velocity of the galaxy expressed in radians/time units; the total angle of 
rotation of the galaxy being a = tmax '■ 



C 



cos {VLq t) sin {VLg t) 
- sin {Qg t) cos {Qg t) 



^2) 



In the astrophysical applications Qg = because the rotation period of a typical 
elliptical galaxy is less than the lifetime of a typical radio source; only in the case of 
NGC326 Qg 7^ 0, see discussion in Section 17^ 

The last translation represents the change of framework from (x*^^)), which is 
co-moving with the host galaxy, to a system (x*^^^) in comparison to which the host 
galaxy is in uniform motion. It should be remembered that the dispersion velocity in 
the cluster (not to be confused with the recession velocity) is ~ 600 km/sec , see 



Table 1 in 



Venkatesan et al. 



()l994f ) . The relative motion of the origin of the 
coordinate system defined by the Cartesian components of the 

galactic velocity vq, and the required 1x3 matrix transformation representing this 
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translation is: 



B 



VJjC 



v-y 



V-Z t 



^3) 



On assuming, for the sake of simplicity, that Vy=0 and Vx=0, the translation matrix 
becomes: 





B 







v_z t 



In other words, the direction of the galaxy motion in the IGM and the direction of 
the jet are perpendicular. From a practical point of view the galaxy velocity can be 
measured by dividing the length of the radio galaxy in a direction perpendicular to 
the initial jet velocity by the lifetime of the radio source; an example of such 
measurement is reported in Section IH^ 

The final 1x3 matrix A representing the "motion law" can be found by composing 
the six matrices already described 



A = B + {C-De-E-F-G) 

(cos {Qg t) sin {Qprec t) sin {"^prec) + siu {Qg t) COS (^prec)) L{t) 

(— sin {Qg t) sin {Qprec t) sin (\l'prec) + COS {Q galaxy t) cos (\E'prec)) L{t) 

V-Z t — sin {"^prec) COS (Qprect) L{t) 



.(85) 



The three components of the previous 1x3 matrix A represent the jet motion along 
the Cartesian coordinates as given by the observer that sees the galaxy moving in 
uniform motion. 
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TABLE 2 

The data from the radio map of NGC1265 



Lrad [ergs/sec] mean spread rate[degree] 



1.110^" 10 



7. APPLICATION OF THE CONTINUOUS MODEL 

The continuous beam model developed in Section IHl is basically represented by 
equations (jKIH) and (jKT|) and can be easily implemented in a code; this is now applied 
to various radio-sources. 

An attempt to reproduce the mechanical power dependence in the jet trajectory is 
carried out in Section 17.51 

7.1. The first part of NGC1265 
The first target chosen for our simulation wa s the radio sour c e NG C1265 at 6 cm. 



Odea fc OwenI (|l986| ) are reported in 



Some basic input parameters extracted from 
Table m The result of the simulation is reported in Figure El and in Figure 0] with the 
radio map superimposed on the simulation ; the adopted input parameters are given 
in the caption. 

The output data as obtained from the code concerning the overall jet length, the 
perpendicular-displacement, and the galaxy velocity of NGC1265 are reported in 
Table El 

7.2. The first part of NGC4O6I 
We have chosen the rad io source NGC4061 at 21 cm, see Figure 5b in 



Venkatesan et al. 



(199J) for a radio map; the basic data as extracted from the radio 
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Fig. 3. Continuous three-dimensional trajectory of NGC1265: the three Eulerian angles 
(English notation) characterising the point of view are 80 ° , Q= 80 ° and 10 °. The 
precession is characterised by the angle *prec= 2 ° and by the angular velocity fiprec= 54.00 
[VlO^year]. The three Eulerian angles are: ° , ■^j= ° and 0^= 90 °. The angle 

of rotation of the galaxy is q;g= °. The physical parameters characterising the jet motion 
are : £^55= 0.02, ti= 10.00 lOVar, no= 0.01 particles/cm^, xi= 18.00 kpc , ri= 0.90 kpc , 
ao= 10 °. 

TABLE 3 

The physical parameters as output of the simulation on NGC1265 



xi [kpc] Vz [km/s] A 2; [kpc] 



18.7 



84.55 



10.3 
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Fig. 4. Figure 1 of lOdea fc Owed (|l98(l ) representing NGC1265 at 4873 MHz superimposed 
on the theoretical trajectory represented by Figure |31 

observations are reported in Table lU The result of the simulation is visualised in 
Figure El and in Figure IHl with the radio map superimposed on the simulation ; the 
caption of the Figure shows the adopted input parameters while the output data are 
reported in Table El 

In this case we know that the dispe rsion velocity in the cluster is 485km/sec, see 



Table 1 in 



Venkatesan et al. 



( 19941 ) . 

The velocity Vz in our model is 119.08km/sec implying an angle of ~ 14.2° between 
the direction of the galaxy motion and the initial direction of the jet. A typical plot 
of the jet velocity as a function of the time is reported in Figure 13 



7.3. The first part of NGC326 



Another test of our code is the radio galaxy NGC326; in particular we refer to the 
image at 1.4GHz as given in Figure 6 (bottom panel) of 



JETS 



27 




Fig. 5. Continuous three-dimensional trajectory of NGC4061: the three Eulerian angles 
(English notation) characterising the point of view are <I>=100 ° , 0=-85 ° and ^=-10 °. 
The precession is characterised by the angle ^prec= 5 ° and by the angular velocity ^}prec= 
60.00 [VlO^year]. The three Eulerian angles are: ° , ^'j= ° and 9^= 90 °. The 

angle of rotation of the galaxy is aG= °. The physical parameters characterising the jet 
motion are : E^q= 0.02, ti= 12.00 lO^year, no= 0.01 particles/cm^, xi= 20.00 kpc , ri= 
1.00 kpc , ao= 10 °. 
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TABLE 4 

The data from the radio map of NGC4061 



Lrad [ergs/sec] mean spread rate[degree] 



0.0610^° 



10 




Fig. 6. Figure 5b of IVenkatesan et alJ ()l994l ) representing NGC4061 at 21cm superimposed 
on the theoretical trajectory represented by Figure 13 



TABLE 5 



The physical parameters output of the simulation on NGC4061 



xi [kpc] Vz [km/s] A 2;[kpc] 



20.9 78.65 14.6 
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Fig. 7. The numerical relationship between velocity and ti in NGC4061. Parameters as in 
Figure [SI 



Murgia. Parma, de Ruiter. Bondi. Ekers. Fanti. fc FomalontI (j200lh : the basic input 



data as deduced from the radio map are reported in Table IHl In this case, the optical 
observations with HST reveal two cores, one coincident with the origin of the radio 
galaxy, and the other at the projected distance of 4.8 kpc. In this case, it is assumed 
that the rotation of the host galaxy is due to gravitational interaction; the total 
angle of rotation being q;g= 90 °. Figure |H1 reports on the obtained simulation and 
Figure IHl reports the radio map superimposed on the simulation ; Table gives the 
data of output from the simulation. 



7.4. A sample of WAT 

A sample of 7 radio-galaxies has already been presented in Section 12.2.11 and is now 
considered as a test of the theory developed in Section El A first look at the 
radio-maps shows a nearly constant opening angle of ~ 5°. On assuming that a 
constant ratio E/xi = -Esg/lkpc = 10~^ characterises the motion, the terminal 
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Fig. 8. Continuous three-dimensional trajectory of NGC326: the three Eulerian angles 
(English notation) characterising the point of view are <I>=110 ° , 0= 65 ° and ^'=-40 °. 
The precession is characterised by the angle ^prec= 2 ° and by the angular velocity ^}prec= 
51.43 [°/10^year]. The three Eulerian angles are: ° , ^'j= ° and 9^= 90 °. The 

angle of rotation of the galaxy is aG=-90 °. The physical parameters characterising the jet 
motion are : E^q= 0.03, ti= 14.00 10'' year, no= 0.01 particles/cm^, xi= 30.00 kpc , rj= 
1.50 kpc , ao= 6 °. 
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The data from the radio map of NGC326 
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Lrad [ergs/sec] mean spread rate[degree] 



0.510^° 



6.45 




Fig. 9. Figure 6 (bottom panel) of iMurgia et al.l (|2f)nih superimposed on the theoretical 
trajectory represented by Figure |H1 



TABLE 7 



The physical parameters output of the simulation on NGC326 



Xi [kpc] Vz [km/s] A 2;[kpc] 



30.6 



98.71 



16.8 
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TABLE 8 

The seven WAT 
ao=5°, rjkpc] = Xi[kpc]/100, 
no = 10^'^particles/cm^ and 
Ese/lkpc = 10-3 



Tannin tkittip 




7?i K rn / 


fi flOVearl 


0647 + 693 


81 


113 


32 


1231 + 674 


35 


263 


6.0 


1333 + 412 


20 


461 


1.95 


1433 + 674 


49 


188 


11 


2236 - 176 


44 


209 


9.5 


3C465 


28 


329 


3.8 


1610-608 


13 


708 


0.82 



velocity Vi and the age of the jets, ti , see Table |S1 can easily be found through an 
iterative procedure: the length of the jet being provided in Table 7 of 



Hardcastle fc Sakelliod (j200J). From Table |H1 the terminal velocities can be 



compared with the central galaxy dispersion velocity in the cluster. Most of the jets 
here considered are distinctly one sided over the length modelled, and this fact is 
here justified by the flip-flop mechanism rather than assuming that they are 
relativistic over the whole length. 
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7.5. Power transition 

The morphological transition between FR-I and FR-II objects can be explained 
under the assumption of constant power P, i.e. negligible synchrotron losses. 
Equation can be modified by inserting the total power P = j and typical 
parameters of the jet previously adopted such as ao = 5° and no= 
0.001 particles/cm^ : 

x(t7)= 36.78 kpc(^^^-^^^y^'' , (86) 

where P33 = P/lO^^Watt. Only the velocity of the host galaxy in the direction 
perpendicular to that of the jet is here considered 

yitj) = 4.05kpc tj f4oo , (87) 

where V400 = f [km/sec] /400. When a comparison of equation (IHUj) and (|H7|) is made, 
it is clear that an increase in the mechanical power of the galaxy corresponds to an 
increase in the elongation along x of the phenomena. A typical plot that tentatively 
reproduces the power transition is reported in Figure El 

There are many 'straight' FR-I's, NGC315 is a good example . This can happen 
when the galaxy's velocity is lower than the standard value here adopted of 400km/s 
or from a particular point of view of the observer. 

8. APPLICATION OF THE TWO-PHASE BEAM 

The two-phase beam developed in Sections HI and El is now applied to two 
radio-sources : in one ,Cygnus A, the relativistic-classic transition is applied and in 
the other ,0647+693, the classic-classic transition. 



8.1. Relativistic Cygnus A 

A good introduction to the cl assical double ra dio sources , as well as to Cygnus A, 
can be found in chapter 6 of Ide Yound ()2002i l . The main physical parameters 
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Fig. 10. Three jet trajectories as a function of the power. P33=100 (dot-dash-dot-dash hne) 
, P33=l (dashed hne ) and ^33=0. 01 (continuous hne). The other parameters in common 
to the three trajectories are ao = 5°, no= 0.001 particles/cm^, Li= 60 , ri= Li/20 , t^=5.0 
and t;4oo=l- The three Eulerian angles are: ° , ° and @j= 90 °. The three 

Eulerian angles (English notation) characterising the point of view are $= 90 ° , 6= 90 ° 
and °. 
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xf^ [kpc] a^^' [°] tfgllWjeai] 



60 5 0.6 



charact erising the jet in Cygnus A , as presented in iPerlev. Dreher. fc Cowan 



( 1984b| ) , are summarised in Table El here xf'^ is the distance between nucleus and 



hot-spots , a^'" is the first opening angle (measured on the radio map) and t°^g , the 
age of the radio-galaxy, as deduced on the basis of the spectral steepening of the 
lobe emission. 

Once the t2 = t°^g^ is inserted in the code, the other input parameters are varied in 
order to obtain this is obtained by solving numerically the integral 

equation (fTSj) . The results of the simulation (as outlined in Section are reported in 
Figure [TT] and in Figure [T21 with the radio map superimposed on the simulation ; 
Table ITUl reports the data. 

8.2. Classic 0647+693 
A radio-map of 0647+693 in the 8.5 GHz band can be found in 



Hardcastle fc Sakellioul (|2004| ) where two radio-images were traced at the resolution 
of 0.3 X 0.6 arcsec and 2 arcsec. The influence of the motion of the host galaxy is 
evident from the maps and can be visualised in the following way 

• A line is traced between the two hot-spots. 

• The perpendicular to the previous line that intersects the nucleus is traced. 
The distance along the perpendicular to the previous line turns out to be 
!^ 36.8 kpc. 
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Fig. 11. Projected 3D trajectory of CygnusA: the three Eulerian angles (English notation) 
characterising the point of view are ^= 90 ° , Q= ° and ^=-20 °. The precession is 
characterised by the angle ^'prec= 1 ° and by the angular velocity ^prec= 600.00 [°/10'' year]. 
The three Eulerian angles are: $j= ° , ^'j = ° and Qj= 90 °. Parameters of the two-phase 
continuous model: xo= 0.0 kpc, ro= 0.02 kpc, /?o= 0.8 , t2= 0.6 10^ year, ti= 0.161 10^ year, 
ao = 5 ° and ai= 35 °. 



JETS 



37 




Fig. 12. Cygnus from the NRAO Image Gallery (white) superimposed on the theoretical 
trajectory represented by Figure ITTI (black) . 



TABLE 10 

The physical parameters output of the simulation on CYGNUS A 



physical quantity\index 





1 


2 


X [kpc] , distance 


0.0 


60.01 


106.18 


r [kpc], rac/iws 


0.02 


5.25 


19.81 


V ^m/s], velocity 


240000 


24577 


6518 


t [lOVear], time 





0.161 


0.6 
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TABLE 11 

The physical parameters output of the simulation on 0647+693 



physical quantity\position 





1 


2 


X [kpc], distance 


0.0001 


80.7 


104.4 


r [kpc], radius 


0.060 


3.6 


8.8 


V [km/ sec] , velocity 


60000.0 


1004.7 


407.6 


t [10'^ year], time 





4.0 


8.0 



From Figure 1 of iHardcastle fc Sakelliod (120041 ) it is also possible to measure the 



second opening angle , ai that turns out to be ^ 25° . 

The results of the simulation (as outlined in Section 0} are reported in Figure IT^ red 
points) and in Figure HU with the radio map superimposed on the simulation ; 
Table ITD reports the data. 

In order to reproduce the long t rails of 0647+693 visible in Figure 1 



of 



Hardcastle fc Sakellioul ()2004( ). the simulation is now performed at a value of t 



which is bigger by factor four in respect to the previous run , see Figure IT^ green 
points) and Figure El where the radio map is superimposed on the simulation. 

9. THE KELVIN-HELMHOLTZ INSTABILITIES 

The very small opening angles of observed extra-galactic radio-jets , if interpreted as 
free expansion, in many cases imply a Mach number M > 6 and existing theories are 
briefly reviewed. Some basic questions can be posed on the stability of such highly 
supersonic jets : 

• Is the jet unstable against perturbation with a wavelength smaller than the 
radius? 
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Fig. 13. Projected 3D trajectory of 0647+693 (red points): the three Eulerian angles 
(Enghsh notation) characterising the point of view are ^= 10 ° , 0= 80 ° and ^'=-10 °. The 
precession is characterised by the angle '^prec= 5 ° and by the angular velocity Qprec= 22.50 
[°/ 10 ''year]. The three Eulerian angles are: ° , ° and @j= 90 °. Parameters of 
the two-phase continuous model: xo= 0.0001 kpc, ro= 0.0600 kpc, c/= 5.0, t2= 8.00 lO^year, 
ti= 4.00 10'' year, oq = 5 ° and ai= 25 °. The green points represents the same trajectory 
after t2= 32.00 lO'^year. 
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Fig. 14. Figure 2 of iHardcastle fc Sakellioul (|200J) superimposed on the theoretical trajec- 
tory represented by Figure IT^ red points). 




Fig. 15. Figure 1 of IHardcastle fc Sakelhoul (120041 ) superimposed on the theoretical trajec- 
tory represented by Figure IT^ green points). 
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• Is the jet stable over long distances? 

• If the jet is unstable, how does the distance over which an instability grows by 
a factor of e depend upon M and the density contrast between the jet and the 
external medium? 

• Are the helical instabilities with an azimuthal number m = 1 and the pinch 
instabilities , m=0 , enough to build a theory on the lateral growth of the jet 
due to an increase of energy in the perturbations ? 

In order to study these questions, we must analyse the Kelvin-Helmholtz instability 
of an axisymmetric flow along the z-axis when the wavelengths A = ^ ( is the 
wave-vector) are smaller or greater than the jet radius a, which is taken to be 
independent of the position along the jet. The velocity , Uq, is assumed to be 
rectangular. The internal ( external ) fluid density is represented by po* (poe) 7 the 
internal sound velocity is s and z/o= The analysis is therefore split in two. 

9.1. Small wavelengths 

If the growth rate of the envelope of the reflecting modes is analysed, it is possible to 
deduce the following approximate relationship for the imaginary {'Jkh) and real {uj.) 
part of the perturbations when M > 1 

s 

7k H = lad- (88) 
a 



with 



7,,^0.17(ln(A;aM))i-2i(^-— )0-06(A;a)0-^ , (89) 



and 



More details can be found in 



uj^^^u^^-^^ka- . (90) 

2 Gi 



Zaninettil (|l986h . 
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9.2. Great wavelengths 



Starting from the equations of motion and continuity, and assuming both fluids to be 
adiabatically cor npressible, it is p ossible to derive and to solve the dispersion relation 



numerically , see 



Zaninetti 



( 19871 ). 

We then start from observable quantities that can be measured on radio-maps such 
as the total length Lobs, the wavelength A^fes of the wiggles (m=l) along the jet, the 
distance A^^^ (m=0) between knots, and the final offset A Lobs of the center of the jet. 
These observable quantities are identified with the following theoretical variables: 



aI, 



Ai 



obs 



\0 ^ \0 



obs 



^oexpl 



2L 



obs 



Mtada 



AL 



obs 



-^obs 



(91) 

(92) 
(93) 
(94) 



where tad = tmin ■ Si/a and Aq is the amplitude of the perturbed energy. The result is 
a theoretical expression for tmin the minimum time scale of the instability, Xmax the 
wavelength connected with the most unstable m ode and /p t he di stance over which 



the most unstable mode grows by a factor e, see 



Zaninetti 



( 19871 ). These parameters 



can then be found through th e set of nonlinear equations prev iously reported. By 



choosing three objects, 



M87 { 



Owen. Hardee, fc Bignelll (1198(1)) . NGC6251 (see 



Perlev. Bridle, fc WillisI ()l984af )l and NGC1265 ( see lOdea fc OwenI (|19860) the 

observational parameters can be measured on the radio image, see Table ^1 

The four nonlinear equations are then solved and the four theoretical parameters are 

found , see Table IT^ 

An application of the results here obtained is reported in Figure El where the 
wavelength of the pinch modes (m = 0) and the oscillations of the helical mode 
(m = 1) as reported in Table IT^ are simulated by identifying the lateral growth of 
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TABLE 12 

Parameters of the observed oscillations in radio-galaxies expressed in averaged 



X Caj\JL± ULiD 


units 






Geoni e f i ■ i rn I ni i s 1 1 ? ' a fi on 


il/87 


A^Gr6251 


A'G'ri265 


Loi,s[averaged radius units] 


48 


58.75 


20 


Xgf^g[averaged radius units] 


4.0 


4.89 


4 


Xlf^g[averaged radius units] 


19.2 


17.5 


13.3 


ALobsidveraged radius units] 


0.96 


5.87 


2 



TABLE 13 

Theoretical parameters from oscillations deduced from the four nonlinear equations 



Theoretical variable M87 NGC6251 NGC1265 



n 


2.75 


2.88 


1.1 


M 


17.28 


19.44 


17.93 




9.61 


28.05 


9.4 




5.57 10-3 


8.2 10-3 


0.25 
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Fig. 16. Parameters as in Figure but with superimposition of the pinch mode (m = 0). 

energy with the precession. Figure IT7I reports the radio map superimposed on the 
simulation of the pinch modes. 

10. SUMMARY 
10.1. New Physics 



The concept of deceleration is both an explanation of observed prope rties of a 



sample of radio-galaxies ,see 



Laing. Parma, de Ruiter. fc Fantil (|l999r i , as well as an 



hypothesis of work to explain the de celeration of a r e. 



ativist ic jet of electrons and 



Laing &: Bridld (|2002bf ) , or by adiabatic models 



positron pairs by mass injection, see 
that requir e coupling between the v ariations of velocity, magnetic field and particle 



density, see 



Laing fc Bridlel (|2002bh . 
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Right AKtniien 1)950.0) 



Fig. 17. Figure 1 of lOdea fc Owed (|l986^ representing NG1265 at 4873 MHz superimposed 
on the pinch mode represented by Figure ITHl 
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Here two new laws of motion for extra-galactic jets in the presence of laminar flow 
are derived , equation (0) and equation (jl3|) . once the two canonical laws of 
resistance are introduced. 

The introduction of the momentum conservation along a solid angle that 
characterises the extra-galactic jet allows us to find: 

• An analytical law , equation (07j), for the velocity and length of the jet as a 
function of the other structural parameters which are: the involved linear 
density of energy, the initial radius, the elapsed time and the solid angle. This 
is the third law of motion. 

• An analytical law for the velocity if conversely the conservation of the 
momentum's flux in a turbulent jet is assumed, equation (jKH|) . This is the 
fourth law of motion. 

• An integral equation (fTHj) of motion if the conservation of the relativistic 
momentum's flux in a turbulent jet is assumed. This is the fifth law of motion. 

The new equation connects the turbulent eddy viscosity with the jet's opening 
angle. 

10.2. Morphological Results 

The introduction of the law of motion allows us to find a space trajectory on 
introducing the jet precession, the velocity and rotation of the galaxy, and the 
inclination of the jet with respect to the plane of the galaxy. 

It can therefore be suggested that t he hydrodynamic al effects necessary t o bend the 



i et (se e for example Section 4.4.5 of Ide Yound (|2002l ) or the discussion in 



Odea 



( 19851 )) can also be explained as an artifact of a change in the system of reference 
combined with precession effects of the jet as well as decreasing jet velocity. From an 
observational point of view it is not easy to discriminate between the bending of the 
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radio-galaxy due to the dynamic pressure , see equation (A35) in Ide Yound (|2002| ) 
and the bending due to the decreasing jet's velocity. Both models require a virial 
speed of the host galaxy ^ 1000^ and differ only in the requirement for the jet's 
velocity which is constant in the traditional approach and variable in the models 
here adopted. This dilemma will be solved when measurements of the jet's velocity 
independent from the radio-models will be available. The spectroscopic observations 
for radio galaxies having ultra steep radio spectra, for example , show velocity shifts 
ran ge from 100 to almost 1000 km/s , 



see 



Roettgering. van Oiik. Milev. Chambers, van Breueel. fc de Kofla (|l997| ). 



10.3. Other approaches 

In FR-I jets, which are thought to decelerate, the profile of this deceleration appears 
in other models that are now briefly reviewed 



In 



Bowman. Leahv. fc Komissarovl ()l996r ) a detailed study of the dynamical 



effects of entraining cool (thermally sub-relativistic) material into hot' 
(thermally relativistic) jets is presented. The dissipation associated with 
entrainment causes only modest loss of kinetic energy flux, and it is shown that 
relativistic jets are affected much less by dissipation than are classical flows. 
Their equation (41) represents the flux of kinetic energy but no analytical 
results are given for the relativistic law of motion such as our formula ()73|). 



In iLaing fc Bridle! ()2002af ) the conservation of particles, energy and momentum 
are applied in order to derive the variations of pressure and density along the 
jets of the low-luminosity radio galaxy 3C 31. Their self-consistent solutions for 
deceleration by injection of thermal matter can be compared with our 
equation (jTfjl . 



In 



Canvin fc LaingI ()2004[ ) some functional forms are introduced to describe the 



geometry, velocity, emissivity and magnetic-field structure of the two 
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low-luminosity radio galaxies B2 0326+39 and B2 1553+24 and an accurate 
comparison between models and data is presented. Their solutions are purely 
numerical and analytical solutions such as our formula (|73|) are absent. 



A comparison 
key-points 



10.4. Analytical results and hydrodynamical codes 
with the hydrodynamical codes is tentatively made on the following 



The hydrodynamical codes are usually performed on a family of parameters 
given by Mf, , which is the beam Mach number , while pb is the inside/outside 
density contrast and t is the time over which the phenomena is followed. Due 
to the fact that the hydrodynamical codes do not give the possibility of closing 
the equations with parameters deduced from the observations, oo figures 
should be produced in order to find the one that be tter appr o ximat es the 
radiosource chosen t o be simulated, see for example iNormanI ()l996l ) and 
Balsara &: NormanI (jl992l ). Conversely, the analytical approach developed here 
, see eqns. fl^. lfT^ . (j^Uj) and (jSH), gives position and velocity once a few input 
parameters are provided. 

The treatment here adopted can not foll ow the nonlinear development of the 



Hardee. Clarke, fc Howelll (|l995|) 



K-H instabilities in 3D , see for example 
where the simulation was performed at Mb=5; but despite this inconvenience, 
the concept of developed turbulence with its consequent constant density along 
the jet is widely used, see Section IHl and Section HI 



10.5. Average density 

Here we have assumed that the density of the jet is constant and equal to its 
surrounding medium. This is certainly something that is not necessarily the case for 
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many astrophysical jets. Some properties can be deduced on the medium 
surrounding the elhptical galaxies from the X-ray surface brightne ss distribution. 



Assuming that the gas is isothermal, the following empirical law 
is used for the electron density: 

^ -1 -3/3/2 



Fabbianol (|l989|)) 



ne{r) = ne(0) 



1+ - 

ax 



(95) 



here /? ~ 0.4 ^ 0.6 , n^{0) ~ 0.1 ^^^^ and ax ~ 1^ 3kpc. On adopting this radial 
dependence for proton density we can evaluate the average density when 
ne(0) = 0.122^^2^, ax = 2.0 kpc and P= 0.4. The average value oi uq , when 
evaluated over 10000 points in the interval 0-100 kpc isn^= 6.74 iQ-SESi^. This 
is the theoretical value of uq that should be used , for example , in equation (jSOI)- 

A. THE ROOTS OF A CUBIC 
To solve the cubic polynomial equation 



aox^ -|- aix'^ + a2X + = 



(96) 



for X, the first step is to apply the transformation 



This reduces the equation to 



X = y — 1/3 — 
ao 



(97) 



+ py + q = 



{9t 



where 



P = 1/3 



Q 2 



q = 1/27 



ao^ 

27 a^ao^ + 2 ai^ — 9 020100 
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The next step is to compute the first derivative of the left hand side of equation 
caUing it /(?/)' 

f{y)' = ^y^ + p ■ (99) 

In the case in which p > , in the range of existence — oo < ?/ < oo , the first 
derivative is always positive and equation has only one root that is real , more 
precisely , 



y = l/6\l-l08q + l2Jl2p^ + 8lq^-2 ^ , (100) 



-108 q+12 y/12 p'^ + 81 q-" 
or equation has the solution in terms of ao,ai,a2 and 03 



1/6. -36^^V^ + 12; 



(3 0200 - ai'^f (3 0200 - ai^f 

4/y h y 



-2/3 (302O0 - oi^) oq'^ ^ .(101) 

(102) 
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